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Abstract. A sample of 229 nearby thick disk stars has been used to investigate the existence of an age-metallicity relation 
(AMR) in the Galactic thick disk. The results indicate that that there is indeed an age-metallicity relation present in the thick 
disk. By dividing the stellar sample into sub-groups, separated by 0. 1 dex in metallicity, we show that the median age decreases 
by about 5-7 Gyr when going from [Fe/H] ss -0.8 to [Fe/H] « -0.1. Combining our results with our newly published a- 
element trends for a local sample of thick disk stars, that show signatures from supernovae type la (SN la), we can here draw 
the conclusion that the time-scale for the peak of the SN la rate is of the order 3-4 Gyr in the thick disk. The tentative evidence 
for a thick disk AMR that we present here also has implications for the thick disk formation scenario; star-formation must have 
been an ongoing process for several billion years. This is further discussed here and appear to strengthen the hypothesis that the 
thick disk originates from a merger event with a companion galaxy that puffed up a pre-existing thin disk. 
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1. Introduction 

Chemical evolution of stellar populations is an important part 
of any model of galaxy formation and evolution. Many studies 
in the past decades show how we are able to further refine and 
constrain models of Galactic chemical evolution by combining 
kinematics and elemental abundances of local dwarf stars (e.g. 
Chiappini et al. l2003l Matteucci l200Tl Edvardsson et al. U993l 
Feltzing & Gustafsson ll998l Feltzing et al. I2003al Bensby et 
al. |2"0"0H Reddy et al. |2UU3}. As evidenced by the cited arti- 
cles our understanding of chemical evolution is mainly driven 
by the studies in the solar neighbourhood but have far reach- 
ing impact for our interpretation of integrated light from other 
galaxies (e.g. Matteucci l2001> . 

However, it is not only the elemental abundances and kine- 
matics of the stars that are of importance when we want to 
further improve the models of galaxy formation and evolution 
but also the ages of the stars (see e.g. Edvardsson et al. 1993 
Raiteri et al. [T996l Pilyugin & Edmunds [T996l Berczik[T333}. 

Many studies have found there to be a clear relation be- 
tween the ages and the metallicities of the solar neighbour- 
hood disk stars (Twarog [T980alll980bl Rocha-Pinto et al. l2000l 
Meusinger et al. 199 lj. In contrast to this Edvardsson et 
al. {1993 1 found no particular evidence for an age-metallicity 
relation in the Galactic disk in the solar neighbourhood and 
Feltzing et al. J200U confirmed this. Feltzing et al. (2001 1 
also showed how dangerous selection effects could be and how 
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an artificial age-metallicity relation can be created (see their 
Figs. 13 and 14). 

Gilmore & Reid ( 1983 1 showed that our galaxy is host to 
two kinematically distinct disk structures. The "new" disk was 
dubbed the thick disk and was found to have a mean metallicity 
around -0.6 dex (Wyse & Gilmore 1995 1 and a scale-height 
of 800-1300pc (e.g. Reyle & Robin I200T1 Chen [T997l while 
the thin disk has a mean metallicity of around -0.1 dex and a 
scale height of 100-300 pc (e.g. Gilmore & Reid [T983l Robin 
et al. 11996b . Recent studies have shown that stars selected to 
belong to either the thin or the thick disk show different trends 
for the elemental abundances (e.g. Fuhrmann 1998 Feltzing et 
al. I2003al Bensby et al. BOOH Bensby et al. 1277011 Reddy et 
al-HHSl Prochaska et al. EOOQI Mashonkina & Gehren lSUOTV 

The question then arises: could it be so that the lack of a 
relation between ages and metallicities for stars in the solar 
neighbourhood is in fact a population effect? That is, are we 
looking at a mixture of stars from (at least) two populations 
that have different chemical enrichment histories? 

It thus appears natural to, yet again, revisit the question of 
an age-metallicity relation in the solar neighbourhood. In the 
study presented here we will address the question of a relation 
between ages and metallicities for stars that are kinematically 
selected to resemble the thick disk closely. 

The paper is organized as follows: in Sect. |5]we describe 
the stellar sample and the kinematical selection criteria and in- 
vestigates if there are potential biases present in the sample. In 
Sect. |3] we describe the choice of ff-enhancement used in the 
isochrones when deriving the stellar ages. In Sect. 0] we derive 
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Table 1. Characteristic velocity dispersions (cry, cr v , and crw) 
in the thin disk, thick disk, and stellar halo, used in Eq. Q. X 
is the observed fraction of stars for the populations in the solar 
neighbourhood and V aS ym is the asymmetric drift (values taken 
from Bensby et al. 2003 and Bensby et al. submitted). 
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Thin disk (D) 


0.90 


35 


20 


16 


-15 


Thick disk (TD) 


0.10 


67 


38 


35 


-46 


Halo (H) 


0.0015 


160 


90 


90 


-220 



ages from stellar isochrone fitting and find that there is a possi- 
ble age-metallicity relation present in the thick disk. In Sect. [6] 
we discuss the implications this tentative age-metallicity rela- 
tion have on the star-formation history of the thick disk, on the 
time-scale of SN la rate in the thick disk, and on our under- 
standing of the origin and evolution of the thick disk. Finally, 
in Sect. we give a short summary. 

2. Stellar sample 

The stellar sample has been selected on purely kinematical 
grounds (see also Bensby et al. 2003, Bensby et al. submit- 
ted). We have assumed that the Galactic space velocities (£/lsr, 
Vlsr, and Wlsr) for the stellar populations in the solar neigh- 
bourhood all can be described by Gaussian distributions. For 
each star (with its specific combination of Ulsr, Vlsr, and 
Wlsr) it is then possible to calculate the probabilities that it 
belongs to either the thin disk (D), thick disk (TD), or the halo 
(H). These can then be used to, for each star, calculate the "rel- 
ative probabilities" TD/D and TD/H. When doing this, the 
fraction (normalization) of the three components in the solar 
neighbourhood must be taken into account. The final relation- 
ship is (see Bensby et al. 2003 1: 
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normalizes the expression; <x w , cry, cr w are the characteristic 
velocity dispersions; V asym is the asymmetric drift; and X is 
the observed fraction of stars in the solar neighbourhood for 
each population. The values for the velocity dispersions and 
the asymmetric drifts are taken from Bensby et al. (2003) and 
the values for the normalizations in the solar neighbourhood 
from Bensby et al. (submitted) (see also discussion below). All 
values are given in Tabled 

A TD/D ratio of 10 means that the star is ten times more 
likely to be a thick disk star than a thin disk star. For our pre- 
vious studies of elemental abundances in the thick and thin 
disks we selected stars with TD/D > 10 as thick disk stars 
and TD/D < 0.1 as thin disk stars. In Bensby et al. (2003 1 we 
also analyzed spectra of a few stars with TD/D between 10 
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Fig. 1. Boxplots showing the distribution in (B-V), My, I/lsr, 
Vlsr, and Wlsr for the full thick disk sample (295 stars, white 
boxes) and the thick disk sample with uvby photometry (229 
stars, gray boxes). In the boxplots the central vertical line rep- 
resents the median value. The lower and upper quartiles are 
represented by the outer edges of the boxes, i.e. the box en- 
closes 50% of the sample. The notches (waists) indicate the 
95 % confidence intervals for the median value. The whiskers 
extend to the farthest data point that lies within 1.5 times the 
inter-quartile distance. Those stars that do not fall within the 
reach of the whiskers are regarded as outliers and are marked 
by solid circles. A common rule is that, samples for which the 
boxplot notch intervals do not overlap are likely to be different 
in their distributions. 



and 2 and it turns out that these stars show the typical thick 
disk trends for the elemental abundances that the stars with 
TD/D > 10 showed. Thus it appears likely that we could use a 
somewhat lower TD/D than 10 in order to increase our sample 
of thick disk stars without compromising the results, but see 
discussion in Sect. [5] 

The most uncertain parameter in the calculation of the 
TD/D ratios is the normalization of the number density of thick 
disk stars in the solar neighbourhood. Different authors quote 
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Fig. 2. The [Mg/Fe] trend for a sample of thick disk stars from 
Bensby et al. ( 2003 1 (circles) and Bensby et al. (submitted) (tri- 
angles). The stars were selected according to the same kinemat- 
ical criteria as discussed in Sect. 13 The solid line indicates the 
degree of a-enhancement that we have used at different [Fe/H] 
(see TableEJ. 



different values: 2 % was found by Gilmore & Reid ( 1983 1 and 
Chen ( flW7t : ~ 6 % was found by Robin et al. JT996I and Buser 
et al. GHUl; and ~ 15 % by Chen et al. ( EoUB and Soubiran et 
al. <U003 ). As there is no simple way to determine which is the 
correct value we inspected colour-magnitude diagrams (CMD) 
of the thin and thick disks, respectively, derived using different 
normalizations for the solar neighbourhood thick disk contri- 
bution (see Bensby et al. submitted). From this inspection it 
is clear that a very low normalization (2 %) must be ruled out 
since then the thin disk CMD showed a population that was 
identical in age and turn-off properties to that of the thick disk 
CMD. The other values of the normalization are, however, not 
as easily distinguished. As a compromise we will use a nor- 
malization of 10 % and a TD/D > 2 and TD/H > 1 as our se- 
lection criteria for thick disk stars. The other parameters, such 
as the velocity dispersions, of course also influence the TD/D 
ratios. They are however better known. 

That the resulting samples are well defined regardless of 
the normalization can be understood through the following 
argument. All stars with TD/D > 2, using a 10% normal- 
ization, will for all normalizations lower than 10% still have 
TD/D > 1. This means that if we select stars with TD/D > 2 
(10% normalization) as our thick disk stars they would also 
be selected as thick disk stars using any of the normalizations 
listed in the literature. In the same way stars with TD/D < 0.6 
(10 % normalization) would always be selected as thin disk ir- 
respective of the value of the thick disk normalization. 

We will restrict the selection of the stellar sample to a cer- 
tain area in the (B - V)-My plane. The limits have been marked 
in Fig. [3] by the dotted lines. This area has been selected on 
the grounds that we do not want to include stars for which 
age determinations are nearly impossible, which is the case the 
lower part of the main sequence where the isochrones are very 
crowded (i.e. My > 6). We also exclude stars that lie on the 
giant branch, since the metallicity calibrations are not valid for 
such stars (see Feltzing et al. 2001 1. In this region the colour of 
the stellar isochrones are also uncertain (see e.g. Yi et al. 1200 II 
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Fig. 3. Example of how a-enhancement changes the stellar 
isochrones. The plotted isochrones are the Y 2 -isochrones from 
Kim et al. (IIUOH and Yi et al. JIqUTI and have ages, [Fe/H], and 
a-enhancements as indicated. The area enclosed by the dotted 
lines are the restrictions in (B - V) and My discussed in Sect.|2| 

their Fig. 2). We have also excluded stars that fall outside the 
following boundaries: (B - V) < 0, (B - V) > 1, and My < 
(see Fig.|3j. 

Selecting all stars in the Hipparcos catalogue that have rela- 
tive errors in their parallaxes less than 25 % (and not flagged as 
binaries or probable binaries), that have radial velocities pub- 
lished in the compilation by Barbier-Brossat et al. ( 1994) we 
get a sample of ~ 12 600 stars. Adding our thick disk criteria 
(TD/D > 2 and TD/H > 1) and the restrictions in (B - V) and 
My given above, we get a sample of 295 stars that are likely to 
belong to the Galactic thick disk. 

Estimates of [Fe/H] for the stars can be derived using 
calibrations of Stromgren photometry (compare Feltzing et 
al. 2001 1. However, not all the stars in the Hipparcos catalogue 
have published Stromgren photometry. 229 out of the 295 thick 
disk stars have Stromgren photometry from the compilation by 
Hauck & Mermilliod (fTWHt . 

We use the metallicities calculated for the catalogue used 
in Feltzing & Holmberg ( EUUPl and Feltzing et al. (EOTTA . In 
Fig. ^ we compare the distributions of (B - V) and My and 
the distributions of the space velocities for the full thick disk 
sample (295 stars) and the thick disk sub-sample which has 
metallicities (229 stars). 

At a quick glance, the distributions of the My, (B - V), and 
the t/ L sR, Vlsr> Wlsr velocities do not seem to be appreciably 
affected when the stars without Stromgren photometry are re- 
jected. We verify this by performing two-sample Kolmogorov- 
Smirnov (KS) tests. The probabilities that the distributions are 
identical is 99.4% for My, 89.7% for (B - V), and >99.9% 
for all three velocities. This means that for all but (B - V) the 
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Fig. 4. CMDs for the thick disk which show how the age distribution increases to older ages as we move to lower [Fe/H]. The 
isochrones are as indicated. We show ages from 20 to 1 Gyr with a step of 1 Gyr. The isochrones for 5, 10, 15, and 20 Gyrs are 
indicated by solid lines. 



distributions are essentially the same. For (B - V) we have in- 
spected the CMDs before and after adding the Stromgren pho- 
tometry. This inspection showed that there is no appreciable 
deterioration in the sampling of red, old stars. This makes us 
comfortable in believing that we do not introduce any of the 
biases discussed in Feltzing et al. (2001 ). 

In summary we do not introduce any bias to our thick disk 
sample by only selecting those stars with published Stromgren 
photometry. 

3. Abundances of n-elements in the stars 

Stars more metal-poor than the Sun have long been known to 
show enhanced levels of a-elements (see e.g. Edvardsson et 
al. U993> . In Bensby et al. (|2003 ) and Bensby et al. (submitted) 
we use the same kinematic definitions of thick disk stars as we 
do here and we derive stellar abundances for a large number of 
thin and thick disk stars spanning a range of metallicities. In 
Fig. |2] we show the resulting [Mg/Fe] vs [Fe/H] trend for the 
thick disk stars in those two papers. 

As can be seen from Fig. [2] the enhancement of [Mg/Fe] 
varies with [Fe/H] for the thick disk. At [Fe/H] = -0.4 the en- 
hancement is ~ 0.3 dex while at [Fe/H] = it has decreased 
to almost solar values. The thick disk stars in Mashonkina 
et al. (2003 1 show the same degree of Mg enhancement, i.e. 
[Mg/Fe] » 0.3-0.4 for metallicities below [Fe/H] = -0.4. (Note 
that the a-element abundance often is defined as the average of 
the Mg, Si, Ca, and Ti abundances, while in Fig.|2]we show 
Mg). 

It is important to know how much enhanced the stars are in 
the ff-elements as a-enhancement has a strong effect on stellar 
evolutionary tracks, and hence on the stellar isochrones that we 
want to use to determine the ages. A set of isochrones with 
an a-enhancement of [a/Fe] = 0.3 is compared to a set of solar- 
scaled isochrones in Fig. [3] The effect of taking a-enhancement 
into account is that ages will be lower than otherwise. 



In the determination of stellar ages from isochrones with 
different [Fe/H] we will use a-enhancements that are in con- 
cordance with the the trend that has been outlined in Fig.|2]and 
is also listed in Table|2] The Yonsei-Yale (Y 2 ) set of isochrones 
(Yi et al. 2001 1 Kim et al. 2002) provides a versatile tool since 
they not only have published calculated sets of isochrones but 
also an interpolator that enables the user to calculate a set of 
isochrones with a specific [Fe/H] and [a/Fe]. 

4. Age-metallicity relation in the thick disk 

Figure|4]shows three CMDs for stars with thick disk kinematics 
and [Fe/H] derived from Stromgren photometry. Each CMD 
is centered at a different metallicity; [Fe/H] = -0.7, -0.4, and 
-0.1, respectively. Stellar isochrones with relevant metallicities 
and a-enhancements (see Sect.|3} are also plotted. 

A visual inspection of these CMDs directly shows that 
the turn-off for the CMD centered at -0. 1 dex is significantly 
younger than the turn-off in the CMD centered at -0.7 dex (see 
Figs.|4k and|4]:, respectively). 

We now proceed to quantify this visual impression. The 
stars are divided into 10 sub-samples according to their metal- 
licities. Each sub-sample has a central metallicity and all stars 
+0.15 dex around this value are included in the sub-sample. 
The central metallicity changes by 0. 1 dex between each sub- 
sample. In this way we create a sliding metallicity binning for 
the age determinations, i.e. the sub-samples are not indepen- 
dent. 

Ages for each sub-sample were simply estimated in the fol- 
lowing way: 

- A set of theoretical isochrones were generated according to 
the central metallicity and appropriate a-enhancement (see 
Sect.fJ. 

- The stars were plotted in the CMD together with the 
isochrones (see examples in Fig.|4}. 

- The plot was inspected and an age was estimated for each 
star (within the restricted area marked out in Fig.fJJi. 
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Table 2. Median ages and spreads for the sub-samples. The 
first column gives the central metallicity, and the second col- 
umn the a-enhancement. Columns 3-6 give the follwoing for 
the samples selected with TD/D > 2: the number of stars 
in the sub-sample, the lower quartile for the age distribution, 
the median age, and the upper quartile for the age distribution. 
Columns 7-10 give the same information but for the samples 
selected with TD/D > 10. 



[Fe/H] [or/Fe] N 



Age (Gyr) 
TD/D > 2 
1/4 1/2 3/4 



N 



Age (Gyr) 
TD/D > 10 
1/4 1/2 3/4 



-0.90 


+0.35 


25 


8.7 


13.6 


15.9 


19 


11.5 


14.0 


16.6 


-0.80 


+0.35 


39 


11.5 


13.5 


16.2 


31 


12.5 


13.5 


16.0 


-0.70 


+0.35 


61 


11.5 


13.2 


16.2 


47 


11.8 


13.6 


16.5 


-0.60 


+0.35 


65 


11.1 


13.0 


15.0 


44 


11.2 


13.0 


14.9 


-0.50 


+0.35 


71 
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12.1 


14.0 


46 
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11.5 


14.0 


-0.40 


+0.30 


58 
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11.7 


13.5 


39 
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13.8 


-0.30 
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-0.20 


+0.15 


60 
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-0.10 


+0.10 


56 


5.5 


7.7 


9.8 


22 


7.0 


8.8 


10.5 


0.00 


+0.03 


35 


4.5 


8.1 


10.1 


14 


6.5 


9.7 


10.5 



- The median age, lower, and upper quartiles where then cal- 
culated for each sub-sample. 

The results are collected in Table [2] The individual stel- 
lar ages are good estimates of the stars age given the set 
of isochrones. More sophisticated methods to derive stellar 
ages from isochrones exist (compare e.g. Ng & Bertelli 1998 
Feltzing et al. 2001 Pont & Eyer 2004 Rosenkilde J0rgensen 
et al. in prep.), However, the age estimates we make are vir- 
tually identical to those from the more sophisticated methods 
(e.g. Rosenkilde J0rgensen who uses a method akin to Pont & 
Eyer 2004 1. Errors on derived ages though are a different matter 
(Rosenkilde J0rgensen private comm.). 

The age estimates together with the central metallicities for 
each bin can now be plotted together in an age-metallicity di- 
agram (see Fig.|5}. We see here clearly that indeed the central 
ages keeps decreasing as the metallicity increases. 

We believe that the change in age with [Fe/H] indeed is 
significant as we are using one set of isochrones and one sam- 
ple of stars that have had their parameters determined in the 
same way. Thus we should not expect systematic errors be- 
tween bins. Furthermore the decline is up to 5 billion years 
from the most metal-poor to the most metal-rich bin. Such a 
large, and systematic, change would be hard to achieve through 
an erroneous analysis. 

5. Thin disk - thick disk confusion? 

In Fig.|5]we see a strong decrease in the mean age for the high- 
est metallicities. Given that the thin disk dominates more and 
more as we go to higher metallicities the question arises: are 
we simply adding more and more of the younger thin disk stars 
into our sample? 

There is in fact some hints in our data that we are picking 
up thin disk stars in our most metal-rich bins. Figure |6] shows 
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Fig. 5. Age-metallicity diagram for the age estimates and cen- 
tral metallicities listed in Table |2]for the samples selected with 
TD/D > 2. We use boxplots to illustrate the age distributions 
within each sub-sample. For a description of box plots see cap- 
tion of Fig. [2 



the distributions of TD/D for the three CMDs in Fig. g] As 
can be seen we are picking up more and more stars with low 
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Fig. 6. Histograms showing the distributions of the TD/D ra- 
tios in three different metallicity ranges for all stars having 
TD/D > 2 (same stars as in Fig.|4j. 
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TD/D ratios as we go to higher metallicities. The stars with 
TD/D < 10 could be intervening thin disk stars. 

Ideally, if the metallicity distribution functions for the two 
disks were well known we could weight the probabilities we 
use to select the stars accordingly. 

The metallicity distributions of the thin and thick disks peak 
at different metallicities (Wyse & Gilmore 11995b . The thick 
disk appears to peak in the interval -0.7 < [Fe/H] < -0.5 
and the thin disk around -0.2 dex. However, the two distri- 
butions are not well constrained. Especially the distribution 
for the thick disk truncates at [Fe/H] « -0.2 in Wyse & 
Gilmore ( 19951. We find that due to the shortcomings in our 
knowledge about the metallicity distribution functions of the 
two disk populations it is not possible to use them to further 
deconvolve the two populations. 

We will here try to address the issue of intervening thin disk 
stars in a different way. In Fig.0we show the same CMDs as in 
Fig-Shut with a stronger constraint on TD/D, i.e. TD/D > 10, 
and compare Fig|6] In this figure we see the same trend (albeit 
with fewer stars) as in the original CMDs, strengthening our 
earlier conclusion that there is an AMR present in the thick 
disk. In Fig. we show the resulting AMR constructed in the 
same way as in Fig. |5] Note, however, in this new case we 
have as a furher precaution excluded all stars with My > 5.4 
in order to make sure that we have as good ages as possible. 
Furthermore, the reader should note that very few stars are on 
the sub-giant branch, where the evolutionary timescale is short 
and hence age determinations can be erroneous (see Pont & 
Eyer 2004). As a further extreme test we have also inspected 
the CMDs where also all star with Vlsr > -50kms _1 have 
been excluded. These CMDs clearly shows that the stars cen- 
tered on [Fe/H] = -0.4 form in the mean a younger population 
that the stars centered around [Fe/H] = -0.7. 

Thus we again find an AMR to be present in the thick disk. 

6. Discussion 

From a spectroscopist's point of view there is now mounting 
evidence that the stars that today have kinematics that we asso- 
ciate with the thick disk show elemental abundance trends that 
are distinct from those seen in stars on solar type orbits (thin 
disk). A recent review can be found in Nissen ( 2003 1. Our own 
results have been published in Feltzing et al ( 2003a 2003b I and 
Bensby et al. ( 20031 |2D04> . see also Fig. [5] These observations 
appear to point to some sort of homogeneity in the chemical 
enrichment process that proceeded the formation of the stars 
that we today associate with the thick disk. 

Our new result that there is a possible age-metallicity re- 
lation in the thick disk, adds a further dimension. The results 
indicate that star formation continued up to 5 billion years in 
the population we identify as the thick disk. 

The observational properties and a standard interpretation 
for the chemical evolution of the thick disk can be described 
in the following way; The abundance trends for the thick disk 
are well-defined and show large a-enhancements at a constant 
value for metallicities below [Fe/H] » -0.4. This is normally 
interpreted as that the star formation was intense and that mas- 
sive stars were the main contributors to the chemical enrich- 
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Age [Gyr] 

Fig. 8. Same as Fig.|5]but using TD/D > 10 instead. 

ment. At higher [Fe/H] the a-enhancement starts to decline to- 
ward solar values. This is typical for what happens when long- 
lived low-mass stars starts to contribute to the chemical enrich- 
ment through the explosion of SN la. The SN la mainly produce 
iron peak elements and none or only little of the a-elements 
which results in a lowering of the [ff/Fe] ratio (compare e.g. 
Fig. 8.6 in Pagei nWfl . 

The time at which the decline starts is not only a function 
of the lifetime of SN la but also of the SN la rate. The exact 
model, and hence lifetime for SN la, is still debated. The two 
most probable scenarios are a double or a single degenerate 
system consisting of either two white dwarfs or a white dwarf 
and a red giant. In the latter case the life time is set by the main 
sequence lifetime of the star that becomes the red giant, while 
in the first case the lifetime for the system prior to become su- 
pernova is set by the time it takes the two white dwarfs to spiral 
in and coalesce, and hence it depends on the initial separation 
and could be longer than the Hubble time. For a recent, in- 
depth discussion of these issues see Livio (2001 1. Furthermore, 
the population synthesis of the SN la progenitors do, in gen- 
eral, not agree too well with the observed frequencies of possi- 
ble progenitor systems (see Livio 2001 and references therein). 
This means that any interpretation of our a-abundance trends 
in terms of a minimum timescale for the star formation period 
in the thick disk will be rather complex (see also Fig. 5.7 in 
Matteucci l2"0"0"TV 

However, our investigation of the relation between ages, 
metallicities, and elemental abundances for stars with typical 
thick disk kinematics points to a conclusion that the SNIa 
rate peaked after a few (~3-4) billion years. We also note that 
the observed abundance trends are tight, compare the oxygen 
trends in Fig. [9] which should indicate that the gas must have 
been rather well mixed throughout. This seems to indicate that 
the gas that the thick disk stars formed out of must have been 
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Fig. 9. Abundance trends for oxygen in the thin and thick disks 
(Bensby et al. 2004 1. Thin and thick disk stars are marked by 
open and filled symbols, respectively. Triangles are data points 
that have been gathered from Nissen et al. (2002i. 

confined to a reasonably "small" physical volume in order for 
the mixing to work efficiently. 

Simulations of galaxy formation indicate that we should 
expect all or a large part of our galaxy to have been assem- 
bled from smaller parts (e.g. Murali et al. 2002 1. The question 
then arises - did these smaller building blocks form stars be- 
fore they merged into today's thick disk? If they did and if 
they were of different sizes, as one may imagine, then they 
would each form their own unique set of elemental abun- 
dance trends which would depend on the initial mass func- 
tion and the star formation rate in each individual building 
block. Then clearly we would not see any homogeneous ele- 
mental abundance trends as has been found by several authors 



(Fuhrmann 1998 Mashonkina & Gehren 2001 Prochaska et 
al. [2000 Feltzing et al. 12003 at . Of course, if these building 
blocks were collected before they formed stars it might be pos- 
sible to find abundance trends in today's thick disk. Another, 
and perhaps more likely possibility, is that what we today ob- 
serve as a thick disk originally was a thin disk that, after having 
produced its stars, was puffed up (see e.g. Robin et al. 119961 
Quinn et al. 1993). Such a scenario would allow the gas to be 



more well mixed at the time of star formation and abundance 
trends would be possible. 

Another possible scenario is a close encounter between the 
Milky Way and another galaxy. This scenario has been mod- 
elled by Kroupa ( 2002 1 who shows that such an event would 
result in a kinematical heating of a pre-existing gaseous disk 
and an increased star formation in this gas. 

7. Summary 

Using a sample of 229 kinematically selected thick disk stars 
we have been able to probe the existence of an age-metallicity 
relation in Galactic thick disk. Ages have been determined in a 
consistent way, using isochrones with the appropriate levels of 
a-enhancement at different metallicities as indicated by recent 
studies that use detailed abundance analysis. Although the stel- 
lar sample is rather small we believe it to be free from severe 
biases that could affect the results. From the investigation we 
are able to draw three main conclusions: 

1 . There is an age-metallicity relation present in the Galactic 
thick disk, indicating that it has had an ongoing star forma- 
tion for a time-period of up to 5 Gyr (this is model depen- 
dent). 

2. The thick disk age-metallicity relation in combination with 
the abundance trends for <?-elements in the thick disk, that 
show signatures from SN la, indicate that the time-scale for 
the peak of the SN la rate in the thick disk is of the order 
3-4 Gyr. 

3. The quite long star formation period in the thick disk 
strengthens the hypothesis that the thick disk formed as a 
result of an ancient merger event between the Milky Way 
and a companion galaxy. 

Studies of the thick disk using nearby stars will always be 
be subject to uncertainties due to the overlapping velocity-, and 
metallicity distributions of the thin and thick disk. At distances 
well above Z* 1.5 kpc from the Galactic plane the thick disk 
is the dominant stellar population. Determining accurate metal- 
licities and (T-abundances for a larger sample of dwarf stars at 
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high Z would therefore enable accurate age determinations that 
would verify the existence (or non-existence) of the thick disk 
AMR deduced from nearby stars. 

Acknowledgements. We would like to thank Poul Erik Nissen and the 
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paper. 

References 

Barbier-Brossat, M., Petit, M, & Figon, P., 1994, A&AS, 108, 603 

Bensby, T., Feltzing, S., & Lundstrom, I, 2003, A&A, 410, 527 

Bensby, T., Feltzing, S., & Lundstrom, I, 2004, A&A, 415, 155 

Berczik, P., 1999, A&A, 348, 371 

Buser, R, Rong, J., & Karaali, S., 1999, A&A, 348, 98 

Chen, B., 1997, ApJ, 491, 181 

Chen, B., Stoughton, C, Smith, J.A., et al., 2001, ApJ, 553, 184 
Chiappini, C, Romano, D., & Matteucci, F, 2003, MNRAS, 339, 63 
Edvardsson, B., Andersen, J., Gustafsson, B., Lambert, D.L., Nissen, 

P.E., & Tomkin, J., 1993, A&A, 275, 101 
Feltzing, S., & Gustafsson, B., 1998, A&AS, 129, 237 
Feltzing, S., & Holmberg, J„ 2000, A&A, 357, 153 
Feltzing, S., Holmberg, J., & Hurley, J.R., 2001, A&A, 377, 911 
Feltzing, S., Bensby, T., & Lundstrom, I., 2003a, A&A, 397, LI 
Feltzing, S., Bensby, T., Gesse, S., & Lundstrom, I., 2003b 
Carnegie Observatories Astrophysics Series, Vol. 4: 
Origin and Evolution of the Elements, ed. A. McWilliam 
and M. Rauch (Pasadena: Carnegie Observatories, 
http : //www . ociw . edu/ ociw/symposia/ series/symposium4 
/proceedings .html) 
Fuhrmann, K., 1998, A&A, 338, 161 
Gilmore, G., & Reid, N., 1983, MNRAS, 202, 1025 
Gilmore, G., Wyse, R.F.G., & Kuijken, K., 1989, ARAA, 27, 555 
Hauck, B., & Mermilliod, M., 1998, A&AS, 129, 431 
Kim, Y.-C, Demarque, P., Yi, S.K., & Alexander, D.R., 2002, ApJS, 
143, 499 

Kroupa, P., 2002, MNRAS, 330, 707 

Livio, M., 2001, Supernovae and gamma-ray bursts: the greatest ex- 
plosions since the Big Bang, (eds) M. Livio, N. Panagia, K. Sahu., 
STScI symp. ser., Vol. 13, p. 334, Cambridge University Press 
lastro-ph/0005344i 
Mashonkina, L., & Gehren, T, 2001, A&A, 376, 232 
Mashonkina, L., Gehren, T., Travaglio, C, & Borkova, T., 2003, 
A&A, 397, 275 

Matteucci, F, 2001, The Chemical Evolution of the Galaxy, 
Astrophysics and Space Science Library, Volume 253, Dordrecht, 
Kluwer Academic Publishers 
Meusinger, H., Stecklum, B., & Reimann, H.-G., 1991, A&A, 245, 57 
Murali, C, Katz, N., Hernquist, L., Weinberg, D.H., & Dave, R., 2002, 
ApJ, 571, 1 

Ng, Y.K., & Bertelli, G., 1998, A&A, 329, 943 
Nissen, P.E., Primas, F, Asplund, M., & Lambert, D.L., 2002, A&A, 
390, 235 

Nissen, P.E., 2003, Carnegie Observatories Astrophysics Series, 
Vol. 4: Origin and Evolution of the Elements, ed. A. 
McWilliam and M. Rauch (Pasadena: Carnegie Observatories, 
http : / /www . ociw . edu/ ociw/symposia/ series/symposium4 
/proceedings .html) 

Pagel, B.E.J., 1997, Nucleosynthesis and chemical evolution of galax- 
ies, Cambridge University Press 

Pilyugin, L.S., & Edmunds, M.G., 1996, A&A, 313, 792 

Pont, F, & Eyer, L., 2004, MNRAS, in press, I astro-ph/0401418 I 



Prochaska, J.X., Naumov, S.O., Carney, B.W., McWilliam, A., & 

Wolfe, A.M., 2000, ApJ, 120, 2513 
Quinn, P.J., Hernquist, L., & Fullagar, D.P, 1993, ApJ, 403, 74 
Raiteri, CM., Villata, M., & Navarro, J.F, 1996, A&A, 315, 105 
Reddy, B.E., Tomkin, J., Lambert, D.L., & Allende Prieto, C, 2003, 

MNRAS, 340, 304 
Reyle, C, & Robin, A.C., 2001, A&A, 373, 886 
Robin, A.C., Haywood, M., Creze, M., Ojha, D.K., & Bienayme, O., 

1996, A&A, 305, 125 
Rocha-Pinto, H.J., Maciel, W.J., Scalo, J., & Flynn, C, 2000, A&A, 

358, 850 

Soubiran, C, Bienayme, O., & Siebert, A., 2003, A&A, 398, 141 
Twarog, B.A., 1980a, ApJS, 44, 1 
Twarog, B.A., 1980b, ApJ, 242, 242 

Yi, S., Demarque, P., Kim, Y.-C, Lee, Y.-W., Ree, C.H., Lejeune, T, 

& Barnes, S., 2001, ApJS, 136, 417 
Wyse, R.F.G., & Gilmore, G., 1995, AJ, 110, 2771 



